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We analyze some aspects of the third law of thermodynamics. We first review both the entropic 
version (N) and the unattainability version (U) and the relation occurring between them. Then, we 
heuristically interpret (N) as a continuity boundary condition for thermodynamics at the boundary 
T = of the thermodynamic domain. On a rigorous mathematical footing, we discuss the third 
law both in Caratheodory's approach and in Gibbs' one. Caratheodory's approach is fundamental 
, in order to understand the nature of the surface T = 0. In fact, in this approach, under suitable 

mathematical conditions, T = Q appears as a leaf of the foliation of the thermodynamic manifold 
' associated with the non-singular integrable Pfafhan form SQrev ■ Being a leaf, it cannot intersect any 

, other leaf S = const, of the foliation. We show that (N) is equivalent to the requirement that T = 

is a leaf. In Gibbs' approach, the peculiar nature of T = appears to be less evident because the 
00 I existence of the entropy is a postulate; nevertheless, it is still possible to conclude that the lowest 

value of the entropy has to belong to the boundary of the convex set where the function is defined. 
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I. INTRODUCTION 



C/5 

O ' We re-analyze the status of the third law of thermodynamics in the framework of a purely thermodynamic formal- 
' ^ ' ism. After a discussion of the status of the third law in current physical literature, and after an heuristic justification 
of the entropic version, we set up a rigorous mathematical apparatus in order to explore the actual necessity for a third 
law of thermodynamics. The approach by means of Pfaffian forms to thermodynamics, introduced by Caratheodory, 
Ph. is the most powerful tool for understanding the problems which can occur in thermodynamic formalism at T = 0. In 
our analysis of the latter topic the Pfaffian form SQrev is expressed in terms of independent extensive variables. One 
j—{ [ finds that T = 0, as an integral manifold of SQrev, can be a leaf of the thermodynamic foliation if sufficient regularity 
J> I conditions for the Pfaffian form are ensured. Contrarily, T = is intersected by the (would-be) leaves S — const. 
C ■ ' which occur at T > 0. The third law appear then as a condition which has to be imposed if a foliation of the whole 
. thermodynamic manifold, including the adiabatic boundary T = 0, has to be obtained. 

Also Gibbs' approach is analyzed. Caratheodory's and Gibbs' approaches together allow to better define the problem 
of the third law. 



The plan of the paper is the following. In sect. II and in sect. Ill a discussion of the third law and of its standard 
proofs is given. A particular attention is devoted to Landsberg's studies, which are under many respects corner-stones 
^ on this topic. In sect. V we remark that Planck's restatement of the third law is not conventional but mandatory 
. for homogeneous systems. In sect. IV, we try to understand, from the physical point of view, if it is possible to 
give a purely thermodynamic justification for the third principle in the entropic version (N). We show that the third 
, principle in the entropic version can be in a natural way interpreted as a continuity boundary condition, in the sense 
P-i that it corresponds to the natural extension of thermodynamics to the states at T = 0. In sect. VI it is shown that, 
^ , in the framework of Caratheodory approach, (N) is equivalent to ensuring that the surface T = is a leaf of the 
thermodynamic foliation associated with the Pfaffian form SQrev ■ The isoentropic surfaces cannot intersect the T — 
surface, because no common point between distinct leaves of the foliation determined by SQrev is allowed. Some 
^ . problems arising when (N) is violated are discussed, and it is recalled that a singular behavior occur if the entropic 
■ ■ ' version (N) fails. In sect. VII a Gibbsian approach to the problem is sketched. We show that the entropy can reach 
its minimum value (if any) only on the boundary T = of its domain. 



*E-mail address: belgiorno@mi.infn.it 
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II. THE THIRD LAW 



The third law of thermodynamics has been formulated in two ways. The original formulation of Nernst concerns 
the behavior of the entropy of every system as the absolute zero of the temperature is approached. Particularly, 
the entropic side of Nernst's theorem (N) states that, for every system, if one considers the entropy as a function of 
the temperature T and of other macroscopic parameters x^,. . . , a;", the entropy difference AtS = S{T, x^, . . . , x") — 
S{T, S\ . . . , S") goes to zero as T ^ 0+ 

lim A5 = (1) 

for any choice of {x^ , . . . ,a;") and of {x^ , . . . ,x"). This means that the limit lim7-^o+ S{T,x^, . . . ,a;") is a constant 
So which does not depend on the macroscopic parameters x^, . . . ,x". Planck's restatement of (N) is 

lim S = (2) 

and it is trivially mandatory for homogeneous systems (cf. sect. V). The other formulation concerns the unattain- 
ability (U) of the absolute zero of the temperature. The (U) side can be expressed as the impossibility to reach the 
absolute zero of the temperature by means of a finite number of thermodynamic processes. Both the above formu- 
lations are due to Nernst, and they are equivalent under suitable hypotheses, as it has been remarked in Refs. [2-4] 
and e.g. also in Refs. [5,6]. 

The third law has a non definitively posed status in standard thermodynamics and a statistical mechanical basis 
for it is still missing. Counter-Examples to (2) have been constructed [7,8], whereas in Ref. [9] models displaying a 
violation of (1) are given. Moreover, the validity of thermodynamics for finite-size systems if T is sufficiently near 
the absolute zero has been questioned. A corner stone of this topic is represented by Planck's objection (see Ref. [10] 
and references therein) against a thermodynamic description of a "standard" system below a given temperature, due 
to a reduction of the effective degrees of freedom making impossible even to define an entropy. The same problem is 
analyzed in Ref. [11] where the breakdown of thermodynamics near the absolute zero is shown in the case of a Debye 
crystal. Thermodynamic formalism is shown to fail because of finite size effects. Indeed, if the finite size of a real 
thermodynamic system is taken into ac;count, according to Ref. [11] near the absolute zero it is no more possible to 
neglect statistical fluctuations in the calculation of thermodynamic quantities like e.g. T, S because they are of the 
same order as the "standard" leading terms^. There is a relative uncertainty in the definition of equilibrium states 
which is of order one. Of course, if one c;onsiders for the number of degrees of freedom a mathematical limit to infinity, 
then the formal success of the thermodynamic approach follows. For more details see Ref. [11]. See also Ref. [12]. We 
don't discuss this topic further on in this paper. 

In Refs. [11,13] it is proposed, in agreement also with the general axiomatic approach of Refs. [2,3], that the third 
law should be assumed as the position of a boundary condition for the thermodynamic differential equations, whose 
experimental validation is stated in regions above the absolute zero. Moreover, according to statements in Refs. [2,4], 
the thermodynamic variables on the "boundary set" of the states at absolute zero temperature could be conven- 
tionally defined as suitable limits (not depending on the path used to approach a particular state at T = 0) of the 
thermodynamic variables in "inner points" of the thermodynamic configuration space and this is proposed as the only 
satisfactory approach to the definition of thermodynamic variables at the absolute zero [2,3]. To some extent, the 
application of the thermodynamic equations to the absolute zero should be considered as a rather formal extrapolation 
of the theory in a region beyond its confirmed domain of validity, and this could be considered as the main reason for 
introducing a new postulate beyond the zeroth, the first and the second law [13]. In sect. IV we come back on this 
topic and give an interpretation of Nernst Heat Theorem as a "continuity" boundary condition for thermodynamics 
at T = 0. 

Concluding this section, it is also remarkable that the third principle, if considered as an impotence principle in 
analogy with the first and the second principle [14], in the (U) version simply does not allow to get T = 0, whereas 
in the (N) version implies also that the work produced by an arbitrarily efficient Carnot machine between T2 > Ti 
(that is, a thermal machine with efficiency arbitrarily near 1"^) vanishes as Ti — > 0^ (see Ref. [15]). For an extensive 
discussion upon the third law see also Refs. [16-18]. 



^The example of Ref. [11] involves a Debye crystal having a volume V ~ Icm^, (N) ~ 10^^; statistical fluctuations are of the 
same order as the leading terms for T ~ 10~*K. 
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III. UNATTAINABILITY VS. ENTROPY BEHAVIOR AT T = 0: LANDSBERG'S ANALYSIS AND 

STANDARD PROOF 



Wc start by discussing (U) and (N) in standard thermodynamics. The double impHcation (U)<^4-(N), according to 
the analysis developed in Refs. [2,4], relies on some hypotheses that it is interesting to recall. 



A. (U)=>(N) in Landsberg's analysis 

A detailed analysis shows that in standard thermodynamics unattainability (U) implies (N) if the following condi- 
tions arc satisfied [2,4]: 

a) The stability condition {dS/dT)^!^ ^.^. > is satisfied for any transformation such that the external parame- 
ters (or deformation coordinates) cc^, . . . , are kept fixed; these transformations be called isometric transformations 
[14]. 

This hypothesis is in general ensured by the suitable convexity/ concavity properties of the thermodynamic potentials 
and is given for ensured in Landsberg's works [2-4] . It is useful to explicit this hypothesis. 

b) There are no multiple branches in thermodynamic configuration space. 

For the condition b) an equivalent statement is "in thermodynamic space no boundary points different from the 
T = ones occur" [2], that is, no first-order phase transitions are allowed. In our setting, this requirement amounts 
to choosing a continuous entropy function. 

c) There is no discontinuity in thermodynamic properties of the system near the absolute zero. 

In Ref. [2] a careful discussion of the conditions to be satisfied in order to ensure (U) is contained. In particu- 
lar, by following Rcf. [2], if a),b),c) hold and moreover (N) fails, then T = is attainable. If a),b) and c) hold, then 
(U) implies (N). If a),b) hold and (N) fails, then (U) implies that a discontinuity near the absolute zero has to occur, 
and such a discontinuity has to prevent the attainability of T = (violation of c)) [2]. Landsberg makes the example 
of an abrupt divergence in the clastic constants of a solid as a conceivable ideal process preventing a solid violating 
(N) to reach a zero temperature state by means of quasi static adiabatic volume variations (the hypothesis of Ref. [2] 
is compatible with the vanishing near T = of the (adiabatic) compressibilities that are related with elastic constants 
in ordinary thermodynamics; particularly, for standard systems one can define the compressibility modulus as the 
inverse of the compressibility; it is proportional to the Young modulus in the case of a solid). Anyway, in standard 
thermodynamics a violation of c) is ruled out and is not discussed further on in Rcf. [2]. Moreover, in standard 
treatment of the third law (U) is associated with the impossibility to get states at T > isoentropic to states at 
T = 0, so that c) is not taken into account. A further discussion is found in the following subsections. 




' > I > 

T T 

(a) (b) 

FIG. 1. (a): Multi-branches structure of the thermodynamic space. According to Landsberg, it implies the validity of (U) 
and the violation of (N). (b): Violation of (N) that implies a violation of (U), due to the presence of the isoentropic AB. 
Landsberg conjectures that (U) holds if a discontinuity near T = occurs. See also the text. In (a) and (b) the dashed regions 
are forbidden. 
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In Rcfs. [2,4] a further condition "entropies don't diverge as T ^ 0+" takes into account the standard behavior 
of thermodynamic systems near the absolute zero. This condition is not necessary if one considers a non-negative 
concave entropy (cf. sect. VII). It can be relaxed when infinite values of the parameters are allowed [19]. E.g., in the 
non-standard case of black hole thermodynamics, the above condition is not necessary, and in Ref. [20] the divergence 
of the entropy occurring in the infinite mass limit for the black hole case is discussed. 

Possible failures of the implication (U)=>(N) are discussed also in Ref. [5], both in the case of reversible processes and 
in the case of irreversible ones. 



B. (N)=>(U) in Landsberg's analysis 

The isocntropic character of the zero temperature states is considered a condition ensuring the unattainability (see 
e.g. Ref. [13] and Refs. [2,3]). A full implication (N)=>(U) is possible in the case of thermodynamic processes which 
consist of an alternate sequence of quasi-static adiabatic transformations and quasi-static isothermal transformations 
(class P(x) according to Refs. [2,3]). Actually, a more general notion of imattainability can be assumed: "zero 
temperature states don't occur in the specification of attainable states of systems" . This is almost literally the (U4) 
principle as in Rcfs. [2,3]. (U4) states that no process allows to reach states at T = 0, even as transient non-equilibrium 
states. Then (N) can fail and (U) can still be valid: In general, the latter hypothesis allows a de-linking of (U) and 
(N) and implies that (N) y^(U) and (U)7^(N) [2,3]. But such a de-linking occurs under particular conditions: the 
failm-e of the implication (U)=>(N) requires again a rejection of one of the hypotheses b),c) above, whereas (N)=>(U) 
fails if processes not belonging to the aforementioned class P(x) allow to reach T = [2,3]. 

For the sake of completeness, we recall in the following subsection also the standard approach to Nernst's theorem, 
which involves heat capacities [13,18]. 

C. (U)-S=>(N) by means of heat capacities: the standeird proof 

The implication (U)=>(N) can be obtained also as follows. It implicitly requires that conditions a),b),c) of Landsberg 
hold. Let us consider two states {Ti,x^ . . . a;"), (T2, y^, . . . , y") and the related entropies 

5(ri,a;\...,x")=5(0,x\...,a;") + j — C,i,...,,.(r) 

S{T2,y\...,yn=S{0,y\...,y^)+J^ _ C,i,...,,n(T) (3) 

where 5(0, x^,. . . , x"), S{0, y^,. . . , y") are the limits of the above entropies as T ^ 0+; in the standard proof one 
assumes also that {Ti,x^ . . . x"), (T2, y^, . . . , y") lie on the same isoentropic surface and that it is possible to perform 
a quasi-static adiabatic process connecting them [13,18]. If T2 = 0, then one gets 

5(0,y\...,y")-5(0,a;\...,x") = y^ _ C,i,...,,„(T) (4) 

and this implies that, if 5(0, y-^, . . . , y") — S{0, x^ , . . . , x") > and the stability condition C.j.i x"{T) > 
holds, a temperature Ti satisfying the last equation always exists, so the unattainability requires S{0, y^, . . . , y") — 
5(0, x^,. . . , x") < 0. The same reasoning applied to the process (T2, y^, . . . , y") — > (0, . . . x") gives the opposite 
inequality 5(0, y^, . . . , y") — 5(0, x^. . . . ,x") > and so one has to conclude that 5(0, y^,. . . , y") = 5(0, x}, . . . , .t") 
[13,18]. The convergence of the above integrals of course requires a suitable behavior for the heat capacities. It 
is remarkable that, according to this standard proof, (U) is implemented by forbidding the presence on the same 
isoentropic surface of states at T = and states at T > 0. Thus isoentropic transformations reaching T = cannot 
exist. This is a key point. Indeed, one could also allow for a different implementation of (U) in which formally 
states at T > and states at T = lie on the same isoentropic surface but, because of some hindrance arising in 
a neighborhood of T = 0, the isoentropic transformation reaching T = actually cannot be performed. This is the 
reason for the hypothesis c) of Landsberg. 

The converse, that is, the implication (N)=>(U), is based on the implicit assumption that processes, which don't 
belong to class P(x) and which allow to reach T = 0, don't exist [2-4]. It is straightforward [18]: Let us consider 
an adiabatic reversible process {Ti,x^ . . . x") — > (T2, y^, . . . , y"). If (N) holds, then there is no possibility to reach a 
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T = state by means of an adiabatic transformation. That is, if (N) holds, then (N) and the second law imply (U). 
Indeed, along an adiabatic transformation 

S{Tux\...,x") = S{T2,y\...,y")^ — C^._,^{T) = _ C,i,...,,n(T) (5) 

and it is evident that, if the final state is (T2 = 0, y^, . . . , y"), then 

^C,i,...,,.(T) = 0, (6) 

which is impossible for C^i > 0. The same conclusion holds if an adiabatic transformation from (Ti, x^, . . . , x") 
to (0, 2/^, ... , 2/") is considered. This proof assumes that the only possibility to get T = is by means of a reversible 
adiabatic transformation. The latter is a reasonable hypothesis, because any thermal contact and any irreversibility 
cannot be successful in obtaining T = due to the second law. For an interesting proof of the above statements see 
also Ref. [21]. 



IV. NAIVE NERNST HEAT THEOREM: A CONTINUITY BOUNDARY CONDITION FOR 

THERMODYNAMICS AT T = 

We assume here a physical attitude, and wonder if it is possible to give a purely thermodynamic justification for 
the third principle in the entropic version (N). This section is dedicated only to an heuristic discussion. A rigorous 
mathematical setting for the third law is found in the following sections. 

We stress that, in our reasoning herein, we adopt substantially Landsberg's point of view as expressed e.g. in Ref. 
[4], p. 69: ". . . one must imagine one is approaching the physical situation at T = with an unprejudiced mind, 
ready to treat a process at T = like any process at T > 0. With this attitude the maximum information concerning 
conditions at T = can be deduced . . ." 

Let us assume that transformations along zero temperature states are allowed. In a reversible transformation at 
T > it is known that {6Q\^^ = T dS. 

As a consequence, 

AS = for adiabatic reversible transformations at T > 0. 

Then, let us consider ideally which behavior is natural to postulate for thermodynamics at T = 0. Along the 
T = isotherm any reversible transformation is adiabatic. From the point of view of thermodynamic formalism, a 
discontinuity with respect to the natural identification between adiabats and isoentropes arises if the states at T = 
are not assumed to be isoentropic. From our point of view, the (N) version of Nernst Heat Theorem appears to be as- 
sociated with a "continuity boundary condition" for thermodynamics at T = 0. Continuity means that the entropy is 
continuous also on the boundary T = and that the identification between isoentropic transformations and adiabatic 
reversible transformations holds also at T = 0. Indeed, by continuity, it is natural, from the point of view of classical 
thermodynamics, to postulate that states at T = are isoentropic and then in the T — S plane the T = line reduces 
to a single point (T = 0, 5 = const). But, in order to match continuously the property that any isothermal reversible 
transformation jt=o at T = is isoentropic, i.e. AS* = V7t=o, one has to require that along the isothermal surfaces 
the entropy variation becomes smaller and smaller, that is, ASt = S{T, x^, . . . , x") — >S'(T, x^, . . . , x") at fixed external 
parameters has to converge to zero as T — »■ 0+. 
The underlying hypotheses are: 

rjo) T = belongs to the equilibrium thermodynamic phase space; 
rji) it is possible ideally to conceive transformations at T = 0; 

772) transformations at T = are adiabatic reversible; 

773) transformations at T = are isoentropic; 

774) there is a continuous match between states at T = and states at T > 0. 

Actually, 774) could even summarize all the hypotheses above, in the sense that a violation of at least one hypothesis 
rjo) ,r]i) ,r]2) and 773) would imply a discontinuity in thermodynamics between zero temperature states and non-zero 
ones. Concerning 771), we recall that Landsberg substantially rejects it, because he postulates a poor population of 
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zero temperature states in order to forbid the T = transformation in the Carnot-Nernst cycle. Each state can be 
associated with its von Neumann entropy and a priori a violation of (N) and a discontinuity are allowed. There is in 
any case a postulate about the density of the zero temperature states which is "discontinuous" with respect to the 
assumptions for the states at T > 0. 

The path of Nernst consists in starting from the violation of the Ostwald's formulation of the second law which is 
implicit in a special Carnot cycle, which has the lower isotherm at T = 0. Sec the figure below. Wc refer to this cyc;le 
as the Carnot-Nernst cycle. If it were possible to perform it, it would imply the existence of a thermal machine with 
efficiency one, which is a violation of the second law of thermodynamics. Note that the violation of the identification 
between adiabats and isoentropes is implicit in the T = isotherm of the Carnot-Nernst cycle. In order to avoid this 
violation, Nernst postulates therefore the unattainability (U) of the absolute zero (see also Ref. [22]). 



A. transformations at T = 

Criticisms against this path, relating the third law to the second one have a long history (see Refs. [3,10,23,24,15,25] 
and references therein) which starts with Einstein's objection. Einstein underlines that near T = dissipations begin 
being non-negligible [23]. This would make the Carnot-Nernst cycle unrealizable because the adiabatic T = could 
not be performed. This kind of criticism could be moved also against any attempt to define transformations at T = 0. 
Nevertheless, it is true that a postulate on thermodynamics is required at T = 0, as variously realized in literature 
(see e.g. Ref. [10]). The objection against the Carnot-Nernst cycle can also avoid referring to irreversibility arising 
near T = 0, as discussed e.g. in Refs. [24,25]. The point is that one reaches the T = surface by means of an adiabatic 
reversible transformation, say BC, and that also any transformation CD at T = has to be adiabatic (see the figure 
below). 




FIG. 2. Carnot-Nernst cycle in the plane T — S. 



Then, it does not seem possible for the system to be carried along the CD transformation contained in the T = 
surface [24,25] because the adiabatic constraint applies to BC as well as to CD and so an operative procedure (no 
matter how ideal) to carry on the cycle seems to be missing. It is to be noted that, by analogy with the operative 
definition of isothermal transformation at T > 0, from a physical point of view, the system should also be considered 
to be in thermal contact only with a "source at T = 0". The non-sense, in the case of T = 0, is evident; a "source 
at T = 0" should be defined (a device able to exchange large (arbitrary) amounts of heat without changing its 
temperature). The point is that any reversible transformation at T = is adiabatic by itself, thus, from the point of 
view of an operative procedure, one should implement an adiabatic insulation of the system. 

This kind of reasoning implies a failure of the thermodynamic formalism at T = 0, because of the impossibility to 
give a satisfactory prescription for implementing transformations at T = 0. In particular, the problem is related with 
the existence of the intersection between adiabatic surfaces (any isocntropic surface intersecting the T = surface 
is an adiabatic surface which intersects the very peculiar adiabatic surface T = 0), because of the apparent absence 
of tools allowing to pass from one to another one adiabatically. In some sense, we find an incompleteness of the 
thermodynamic formalism at T = 0, because there are serious problems in defining an operative procedure [23-25]. 
Nevertheless, we wish to underline that, in line of principle, it could be still possible to implement the adiabatic 
transformation at T = as a distinct adiabatic transformation , because, even if an adiabatic constraint is required. 
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it corresponds to a path mathematically distinguished in the thermodynamic space. It is clear that, if one considers 
a system described by {U, V, N) and, in ideally approaching T = considers the system as closed, then no possibility 
to distinguish operativcly between the adiabatic transformation implemented in order to approaching T = and the 
adiabatic transformation T = is left. But for systems with a larger thermodynamic space (e.g. systems characterized 
by other deformation parameters) one could have closed systems where the adiabatic and isothermal transformation 
at T could be implemented. 

A final comment. It is often stated that the third principle is not as fundamental as the first and the second 
ones and that it is not related to any new potential in thermodynamics, whereas the first law is associated with the 
internal energy function and the second law with the entropy [5]. On this side of the topic, the third law prescribes 
the behavior of the zero-temperature part of the entropy limT^o+ S{T, x^, . . . ,x'^) = S{Q, x^, . . . ,x") = , . . . , a;") 
and fixes its value to zero. We shall show that 

The third law of thermodynamics (Planck restatement) corresponds to a regularity condition of the Pfaffian equation 
5Qrcv = on the boundary T = of the thermodynamic manifold. It is equivalent to the request that a well-defined 
foliation of the whole thermodynamic manifold exists. 

In the following sections, we discuss the problem in a mathematically rigorous framework. 



V. ABSOLUTE ENTROPY AND PLANCK'S POSTULATE 



In our discussion of the third law, the zero-tcmpcrature entropy constant is undetermined, with the only constraint 
Sq > suggested by statistical mechanical considerations. Planck's restatement of (N) requires 5*0 = 0, that is, 
S' — > 0+ as T ^ 0+, because the constant So (entropy at T = 0), which does not depend on the thermodynamic 
parameters, does not affect any physical measurement [26]. According to some authoritative experts in the field 
of thermodynamics, this corresponds to a sufficient condition for implementing (N), not a necessary one [17,15,27]. 
Problems arising with chemical reactions can be suitably solved [17,15,27]. An analogous position against the neces- 
sity of Planck's restatement appears also in statistical mechanics. E.g. in Ref. [28] statements, according to which 
systems violating S* ^ 0+ as T ^ 0"*" in the thermodynamic limit, automatically violate (N), have been criticized. 
Residual entropies coming from theoretical calculations in statistical mechanics, as far as they are not involved with 
a dependence of the ground state entropy on macroscopic parameters, they still cannot be considered as violations 
of the third law [28]. Also in Ref. [29] the violation of Planck's statement is not considered a priori as implying a 
violation of (N) . We first discuss the problem in the framework of the thermodynamics of homogeneous systems; then 
we add some comments about the relation with statistical mechanics. 



Notice that a necessary condition for (N) to hold is that S is continuous in the limit T = 0, whichever state is 
considered on the surface T = 0. In fact, let us define X = x^,. . . ,x" and let us assume that S is not continuous 
in (0,Xo). If this discontinuity is not simply an eliminable one, there exist two different sequences {Tn \Xn^}, with 
i = 1,2, such that {Tn\Xn^) (0, Xq) as n — » oo and, moreover, such that 

lim 5(r(i),X(i)) ^ hm 5(rf ),X(2)). (7) 

(N) is badly violated. See also sect. VII. The violation of (N) occurs also in the case of a (unnatural) eliminable 
discontinuity. As a consequence, the continuity of S on the surface T = is assumed. 

It is to be noted that, if (N) holds, the entropy constant at T = cannot depend on the composition variables 
n* which specify the number of moles of the component substances which are present in the material whose thermo- 
dynamic properties are studied. Herein, we let composition variables to be included in the set of what we called 
deformation parameters ^. Then, under suitable hypotheses it holds 

lim -^=0. (8) 

T^o+ an' 



^This choice can be opinable in light of a rigorous axiomatic approach [30], but it allows us to call deformation parameters 
all the parameters different from U (from T) in our discussion, which is limited to some aspects of the third law. 
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This can be deduced also by means of homogeneity properties of the entropy; for a pure phase at constant p, T one 

has 



dn 

Moreover, one has for a fc-components system at constant p, T 

k OCT 

SiT,V,...,n\...,n'^) = J2n'^ 



S = n'' (9) 



9n' 



i=l 

k 



(Note that is not the entropy of the single component; such an identification would originate a wrong expression 
for the entropy, as it is clear for the case of mixtures of ideal gases, where an entropy of mixing appears). 
As a consequence, it appears that the arbitrary constant So is zero both for pure phases and for mixtures and chemical 
reactions. In the latter case, the third law states that 

lim Y ^=0 (11) 

i=l 

where z/* are the stoichiometric coefficients. Actually, each derivative should vanish at the absolute zero. The 
alternative definition of Ref. [27,15] seems to be not satisfactory from the point of view of (N), because the zero- 
temperature entropy appears to depend on composition variables (which would allow for the composite a different 
zero-temperature entropy for different molar fractions of the components, against the postulate of Ncrnst). 
Then it is thermodynamically appropriate to put 5o = 0. This actually not only does not amount to a real loss of 
generality, because measurements leave the constant undetermined it is also a necessity (if the third law holds) for 
homogeneous systems, because So is required to satisfy the homogeneity property of the entropy of a homogeneous 
system. S is an homogeneous function of degree one in the extensive variables, say, {U, V, N). Because of the Euler 
theorem, this implies that, by introducing the operator 



the entropy satisfies the equation 



by inverting 



Y S = S. (13) 



dS_^ 
dU 



= g{U,V,N) (14) 



with respect to U [which is allowed by the fact that dT/dU > 0] one finds U = h{T,V, N) which is an homogeneous 
function of degree one in the extensive variables {V, N). Then one obtains S{T, V, N) which is homogeneous of degree 
one in the extensive variables {V, N) [it is a quasi-homogeneous function of degree one and weights (0, 1,1)]. If (N) is 
satisfied, then 

lim S(T,V,N) = So (15) 
for any choice of V, N. Because of the homogeneity, for any A > one has 



^In order to understand this point, it is important to underline that the constant So has actually no operative meaning, in 
the sense that thermodynamic measurements (and extrapolations for the limit T — > 0"*") are relative to the integral of C/T. 
So, in line of principle, it can be put equal to without affecting thermodynamic measurements. 
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lim 5(T, AV, AA^) = A lim S{T,V,N) = XSq 



(16) 



which is consistent with the independence of the Hmit from V, N only for = 0. No additive constant can appear as 
T ^ 0+, because of the homogeneity, thus Planck's restatement of the third law is mandatory if (N) is satisfied. 

Summarizing: 

Planck's restatement is mandatory if (N) holds, due to the homogeneity property of S. 

A short comment about the third law in statistical mechanics is in order. It is commonly stated that a violation 
of (N) occurs if the ground state is degenerate. Moreover, As far as the limit as T — » 0+ is concerned, one has to 
distinguish between finite size systems and bulk systems. In the latter case, Griffiths shows that, in determining 
the behavior of thermodynamic systems near the absolute zero, in measurements, what is really important is the 
contribution of the excitable low-energy quantum states: for bulk systems the contribution of the ground state, at 
reachable low temperatures, is irrelevant in determining the behavior of the system, which is instead dominated by 
the contributions of the low-lying energy states [7] (the degeneracy of the ground state is not a good indicator of the 
behavior of the entropy for bulk systems at low temperature because the thermodynamic limit has to be carried out 
before the limit as T — > 0+ and the two limits in general don't commute [7]). In statistical mechanics, the ground state 
degeneracy for bulk systems does not play a straightforward role in determining the behavior as the absolute zero is 
approached, and examples exist where the ground state is not degenerate but the limit S* ^ 0"*" is not implemented [7] . 
However, a role for the degeneracy of the ground state can be suitably resorted as in Ref. [8]. Therein, it is remarked 
that the entropy functional at T = depends on the boundary conditions. Different boundary conditions correspond 
to different ground states for the bulk system, and the contribution of the excitations of the low-lying states near the 
absolute zero can be related with a maximally degenerate ground state by means of a variational criterion [8] . 
We limit ourselves to refer the reader to Ref. [31] for a further approach to the problem of the third law by means of 
the concept of dynamical entropy and to Ref. [32] for another interesting point of view concerning the problem of the 
third law in presence of ground state degeneracy. 



In Caratheodory's approach [33], the infinitesimal heat exchanged reversibly SQrev^ defined on a open simply 
connected domain V, is a Pfaffian form, i.e. a one-form w, whose integrability has to be ensured in order to define an 
entropy function. See e.g. Refs. [14,34 37]. This approach appears to be very clarifying with respect to the problem 
represented by the special surface T = 0. In the following, we use w = SQrev 



Caratheodory's principle of adiabatic inaccessibility is usually stated for the case where V has no boundary, that 
is, dX) = [36,37]. It can be formulated as follows: 

(C): eacii neighborhood of any state xq belonging to the domain V contains states which are inaccessible from 

xq along solutions of lo = 0. 

This principles ensures that the Pfaffian form w is completely integrable, i.e. it satisfies u A du = 0, in such a 
way that a foliation of the thermodynamic manifold into isoentropic hypersurfaces is allowed. 
If a boundary is present, there are some changes in the theory^. The integrability condition 



has to be imposed in the interior of the domain of the differential form w, where oj is required to be at least C^. These 
properties ensure that Probenius theorem can be applied and one obtains a foliation in the inner part of the manifold. 



"^The author is indebted to Lawrence Conlon for an enlightening e-mail about the problem of Frobenius theorem for manifolds 
with boundary. 



VI. CARATHEODORY'S APPROACH AND (N) 



A. Foliation in thermodynamics 



u> Ad(j = 



(17) 
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For what concerns the boundary, it can be in part transverse and in part tangent to the inner foliation. It is tangent 
when it is a leaf of the foliation itself, i.e. if the boundary is an integral manifold for w [38]. If, instead, it is not a 
leaf, one can induce on the boundary a foliation from the inner foliation. Then, a foliation of the whole manifold is 
obtained if sufRcient regularity conditions for u) on the boundary are assumed. 

Let us now consider what happens in thermodynamics. The integrating factor T vanishes at T = 0, which means what 
follows. The non-singular integrable Pfaffian form SQrev gives rise to a foliation of the thermodynamic manifold for 
T > 0. Each leaf of the foliation is a solution of the equation dQrev = 0. This foliation has codimension one (i.e., each 
leaf is an hypersurface in the thermodynamic manifold). For T > 0, the leaves of the foliation are the hypersurfaces 
S = const. One has then to determine if the surface T = is a leaf itself. It is indeed an integral submanifold of 
the Pfaffian form lo, in the sense that any curve contained in the surface T = 0. For any initial point lying on the 
submanifold T = 0, there is a curve 7 lying entirely in the submanifold T = 0. One has to ensure the uniqueness of 
the solutions of the Cauchy problem for the ordinary differential equations associated with lu. The Lipschitz condition 
for each of them would be enough in order to get an unique solution. If uj is also on the boundary T = 0, then 
we can show that the uniqueness is ensured and T = is a leaf of the thermodynamic foliation (a tangent leaf). The 
special leaf T = cannot intersect any other leaf S = const, defined at T > 0, because no intersection of leaves is 
allowed. In the following subsections, we analyze the above problem in detail. 



B. domain D 



Let V be the thermodynamic manifold whose independent coordinates are the extensive variables U,V,X^, . . . , X"; 
the variables V, X^, . . . , X" will also be called deformation parameters. Assume that diml> = n+2. V is assumed to be 
a open convex set, in order to match the concavity property of S. Homogeneity requires that (A U, X V, X , . . . , A X") 
belongs to V for each real positive A, thus V has to be also closed with respect to multiplication by a positive real 
scalar, i.e., V has to be a cone. Then, it is natural to require that "D is a convex cone [39]. One can also relax to some 
extent the latter condition [e.g., a positive lower bound on V, A'^ should be introduced on a physical groimd, they cannot 
be arbitrarily near the zero value or statistical fluctuations would not allow to define a meaningful thermodynamic 
state. Of. [39]]. 



C. Pfaffian forms and homogeneous systems 



Let oj = 6Qrev be the Pfaffian form of interest, which is identified with the infinitesimal heat exchanged reversibly. 
It is assumed to be at least of class in the inner part of the thermodynamic manifold. One can write 

SQrev =dU+pdV -Y, ^idX\ (18) 

i 

where {U, V, X^, . . . , X") are extensive variables. The integrability of SQrev ensures that 

SQrev = T dS. (19) 
We assume that SQrev is an homogeneous Pfaffian form of degree one. This means that the vector field 

i 

is a symmetry for SQrev [40,41], in the sense that 

Ly SQrev = SQrev: (21) 

where Ly is the associated Lie derivative. It can be shown that, in the homogeneous case [39], an integrating factor 
for (18) exists and it is given by 

/ = iySQrev = SQrev{Y) = U + p V - ■ (22) 

i 

The integrating factor is required to be such that / ^ 0, which means that Y is not a characteristic or trivial symmetry 
for the distribution associated with SQrev Cf. Ref. [41]. Moreover, one requires / > 0, which is easily shown to be 
equivalent to the conventional choice T > 0. We sketch here some results of Ref. [39]. One finds 
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SQrev = fdS (23) 

and it can be shown in general that, for any homogeneous integrable Pfaffian form, cj/f has to be equal to dH/H, 
where /J is a positive definite homogeneous function of degree one; moreover, the homogeneous function H is unique 
apart from a multiplicative undetermined constant [39]. This function H is actually the entropy S, as it can be 
straightforwardly deduced also by direct comparison with the definition of S as extensive function 

S=^U+^V-J2 (24) 

i 

in fact, one finds that / coincides with the product TS. As a consequence, one has 



which implies 



dS=^ = % (25) 



So exp( j), (26) 



where F indicates a path between a reference state Uq, Vq, Xq, . . . , Xq and the state U, V, X^, . . . , X". We require that 

the thermodynamic foliation is described everywhere in by the leaves S = const., which means that S has to be 
defined everywhere on the thermodynamic manifold (except maybe on the boundary / = 0) [39]. The only problems 
can occur where / = 0. Moreover, one also assumes that to each level set S = const, correspond a unique leaf (which 
means that each isoentropic surface is path-connected, as it is natural to assume). 

D. zeroes of the integrating factor and the domain 

Let us define the set 

Zif) = {iU,V,X\...,X^)\fiU,V,X\...,X") = 0}. (27) 

Z{f) is the set of the zeroes of /. Wc define also 

Z{T) = {{U,V,X\...,X^)\T{U,V,X\...,X")=0}, (28) 

and 

Z{S) = {{U,V,X\...,X^)\S{U,V,X\...,X") = 0}. (29) 
The set Z{f) = Z{T) U Z{S) corresponds to an integral manifold of w (w is non-singular and w = fdiS). 

1. Z{T) 

The set Z{T) is expected to be an hypersurface, but, in general, it could be a priori a submanifold of dimension 
1 < A; < n -|- 1. Actually, it is natural to assume that it is a hypersurface, i.e., a manifold of codimension one. The 
equation 

T([/,V,X\...,X") = (30) 

is required to be implemented for any value of V, X^, . . . , X" which is compatible with the system at hand. Contrarily, 

one should admit that T = could be allowed only for a restricted region of parameters (e.g., a crystal could not 
be allowed to assume a value V — Vq for the volume at T = 0) in such a way that a thermal contact with a lower 
temperature system could not lower the system temperature near the absolute zero if values of the parameters outside 
the allowable range would be involved. We then assume that the T = is a path-connected hypersurface which 
coincides with the adiabatic boundary T = of the thermodynamic manifold. See also subsect. VI F . 
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2. Z(S) 



The set Z{S) has to be contained in the boundary of the thermodynamic manifold. This is a consequence of the 
concavity of 5* and of the requirement S > 0. as it is shown in sect. VII. 5 = can be moreover attained only on 
the boundary surface T = 0, in fact 5 = at T > can be rejected on physical grounds. In fact, any state z such 
that Tz > and S{z) = should have the peculiar property to allow the system only to absorb heat along any path 
7^ starting from z in a neighborhood of z. If Cj{T) is the heat capacity along a path 7 which does not contain 
isothermal sub-paths, one has that 



(T) 



(31) 



should be positive for any state y non isoentropic to z in Wz, which is possible only for heat absorption (in fact, 
C-y^ (T) < would be allowed for states such that <Ty, which would imply heat absorption, and Cj^ (T) > would 

be allowed for states such that > Ty). The same is true if one considers an isothermal path starting at z, in fact 
the heat exchanged would be T^AS and AS* should be positive in a neighborhood of z, being S = a global minimum 
of S. Thermal contact with a colder body at T < should allow an heat flow outgoing from the system because of 
the second law in Clausius formulation. Then, no possibility to approximate such a thermal contact by means of a 
reversible transformation exists, and this behavior can be refused as pathological. 

There is also another argument one can introduce against the possibility that, for a non-negative definite entropy, the 
set Z{S) — Z{T) is non-empty. By using standard formulas of thermodynamics, one has 



S{T,x\...,x"+') = S{0,x\...,x"+')+ _C,i,...,,.+i(z) >0 VT; 

Jo z 



(32) 



{x^, . . . ,x"^^) are deformation parameters (they could be also in part intensive); because of the concavity condition 
Cx^,...,x"+^{T) > the entropy can vanish only for T = 0. A further discussion is found in subsection VI F. 

Then a non-negative concave entropy implies that the set Z{S) of the zeroes of S is contained in the set Z{T) of 
the zeroes of T: 



Z{S) C Z{T). 



(33) 



The two sets coincide if (N) holds, otherwise Z{S) C Z{T) and it could be that Z{S) = 0. Then we get the following 
equality: 



Z{f) = Z{T). 



(34) 



If one considers a concave entropy which can be also negative, then it happens that Z{f) D Z{T) because Z{S) is 
not, in general, a subset of Z{T). A typical example is the classical ideal gas. Let us consider the monoatomic ideal 
gas. One has [48] 



S{U, V,N) = N 
the corresponding Pfafiian form is 



5 f U^^V 1 
2 + V N^^^ (3 7r)3/2 



2 U 2 [/ , /C/3/2 V 1 



dN 



(35) 



(36) 



and one has 



and 



•' 3 



T = 



2 U 

3 TV' 



5 / U^I^V 1 

2 V ^^/^ (3 7r)3/2 



(37) 



(38) 



In this case one has 
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Z{T) = {[/ = 0} (39) 



and 

Z{S) = {{U,V,N) I = exp [--p. (40) 

Then Z{f) D Z{T) and / vanishes before [/ = is reached. 

The equation / = is an impUcit equation which defines a submanifold of the thermodynamic manifold. This is 

trivial if / is at least everywhere in VUdV, in fact / = defines a hypersurfacc contained in the domain. This 
submanifold could be trivially an hyperplane U = Uq = const., or a non-trivial hypersurface U = b{U, V, X^, . . . , X"). 
A further discussion is found in subsect. VI F and in subsect. VII C. 



E. T = in thermodynamics 

The surface T = is usually excluded as unphysical. In Rcfs. [30,42,43] an axiomatic approach to thermodynamics 
excludes a priori that the value T = belongs to the thermodynamic manifold. In literature often the claim appears 
that the temperature, for consistency, has to be strictly positive (see e.g. Refs. [42,43]). From the point of view of 
the approach involving extensive variables, it has to be discussed if Z{T) is empty or not (it is surely non-empty in 
the black hole case, sec also Rcf. [44]). In the former case, it should be discussed if there is a lowest temperature 
[2] (maybe different for each system) and what this implies for the physics. A different lowest temperature for each 
system is not a viable hypothesis, because one could put in thermal contact a system at its own lowest temperature 
with another system at a lower temperature and should see an heat flux from the former to the latter, and a decrease 
of the temperature of the former. Thus, a lowest temperature should be allowed to be the same for all systems. 
Moreover, from an experimental point of view there is no apparent limit to the possibility to approach T = 0. From 
a theoretical point of view, there is actually no physical hindrance to consider T = as a possible value. From a 
mathematical point of view, the Pfaffian form {SQrev) vanishes at T = but it is non-singular. A singularity of a 
PfafRan form uj = ai{x) dx^ is defined as the set where ai{x) = i, i.e., where all the coefHcients of the PfafRan 
form vanish. But in the thermodynamic case, no singularity is allowed, because of the coefficient of the internal energy 
term, which is in any case one. Thus, no mathematical hindrance to consider T = in the thermodynamic domain 
appears. The singularity appears only when T is used as independent variable, and it is due to the fact that the 
change of variable U ^ T \s a diffeomorphism only for T > 0. See also the following subsection. This topic is also 
discussed in Ref. [19]. 

F. boundary revisited 

In thermodynamics, as discussed in sect. VI D, it is to some extent natural to assume that the boundary T = is 
described explicitly by a (maybe smooth, let us assume at least C^) function: 

C/ = 6(X\...,X"+i); (41) 

one can figure that it corresponds to the equation for the ground-state energy of the system as a function of the 
deformation parameters, as it is clear from the following analysis. 6 is a function which is homogeneous of degree one 
with respect to (X\ . . . , X"+i): 

6(AXS...,AX"+i) = A6(X\...,X"+i). (42) 
Thus, 6 has to be defined on a cone K,}, C IR""*"^. Moreover, if Uq, X^, . . . , belongs to the boundary T = 0, from 

dU —{U,X\...,X"+^), (43) 

where the integral is an improper integral, because dT/dU = l/Cxi,...,xn+i — > oo as T ^ 0+, and from the concavity 
of S, which implies that Cxi,....x"+i > 0, one finds that U > Uq, i.e. it has to hold U > b{X^, . . . ,X"+i). Thus, 
the domain V has to be such that the inequality U > b{X^ , . . . , X'^'^^) is implemented for each U and for each 
(X\ . . . , X"+i) e tCb. The domain V contains the set 
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epi(6) = {{U,V,X\...,X^)\{V,X\..., &K^,U> h{V, X\..., X")}. (44) 

This set is the so-called epigraph of the function b. If the function b{V,X^, . . . is required to be convex, then 

it is defined on the convex cone /Cb, and its epigraph epi(6) is a convex cone (the epigraph of an homogeneous 6 is a 
cone). Then, the domain V can be chosen to be 

I? = epi(6). (45) 

One can also assume that 2? is a convex cone of the form 

V = {{U,V,X\...,X^) I (y,X\...,X"+i) G/y x7^i x...Ixn, U>b{V,X\...,X^)}, (46) 

where the intervals Iv,Ix^ j • • • j ^X" are IR-|-. Wo can find a coordinatization of the boundary by means of coordinates 
{B, X^, . . . , X"+^) such that the boundary T = coincides with B = 0. In fact, we can simply define 

B = [/-6(X\...,X"+i); (47) 

B > is a degree one homogeneous function, and dU/dB = 1. By inverting one finds 

[/ = B + 6(X\...,X"+i). (48) 

As a consequence, one gets 

f = f{B,X\...,X^+')=B + b{X\...,X^+^)-Y, ^kX"; (49) 

k 

by definition, / vanishes for B = 0, i.e. 

= b{X\ X"+i) - ^ a(0, . . . , X"+i) X'^ (50) 

k 

^=a(0,X\...,X"+l) (51) 

Notice that, by defining for all i = 1, . . . , n + 1 

f)TT 

^i{B,X\..., X"+i) =^i{B,X\..., X"+i) -—{B,X\..., X"+i), (52) 

one finds 

(jj^dB-J2 ii{B,X\...,X''+^) dX\ (53) 

i 

and it holds ^i{B = 0,X^, . . . ,X"+^) = for all i = 1, . . . ,n + 1, because of the definition for / = to be an 

integral hypcrsurfacc for lu. Moreover, notice that, under this assumption about the boundary T = 0, one obtain that 
Z{S) C Z{T) necessarily. In fact, one can write for an everywhere continuous entropy 

s{B,x\...,x-+^) = s{o,x\...,x-+^) + dy T(Y,x\...,x-+^y ^^^^ 

where S'(0, X^, . . . , X^^^) is the value attained by S* at _B = by continuity; it is evident that S cannot vanish outside 
Z{T), because 5(0, X^, . . . , X"+i) > and /^^ dY l/T{Y,X^ , . . . ,X''+^) > for all B > 0. 

As far as the entropy 5 as a function of B, X^, . . . , is concerned, it is such that 

dS 1 



dB T(B,Xi,...,X"+i)' 
dS UB,X\...,X^+^) 



dxi r(s,xi,...,X"+i) 



(55) 

Vi = l,...,n+l. (56) 
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One could also allow for different choices of coordinates. Instead of B defined as above, one could introduce another 
(maybe local) coordinate for the boundary / = 0, say B, such that the boundary coincides with B = 0. The coordinate 
transformation U B is required to be regular, i.e. 

dU 

-^j^O for S = (57) 
dB 

BB 

— ^0 iovB = 0, (58) 



then one can find 



where 



r)TI 

u; = ^dB-Y,UB,X\..., X"+i) dX\ (59) 



and 



ii{B,X\ X"+i) ^ ^i{B, X\ . . . ,X"+i) - ^{B,X\ X"+i), (60) 

/=^^-E^^(^'^''---'^"^')^'- (61) 
oB 

It is important to point out that, also in these coordinates, one has 

|i(^ = 0,X\...,X"+i) = 0. (g2) 

In the following, we define 

g^a(B,Xi,...,X"+i). (63) 
oB 

Notice that, because of the properties of the Pfafhan form cj, the absolute temperature T cannot be used as a good 
coordinate for the boundary, in fact dU/dT — > as T — > 0+ for all physical systems allowing a finite 5 at T = 0. This 
choice (as well as the choice of /) seems to transform the regular PfafRan form w into a singular one, but this trouble 
is simply due to the singularity in the jacobian of the coordinate transformation U >—> T, which is a diffeomorphism 
only for T > 0. See also [19]. 



G. condition to be satisfied in order that T = is a leaf 



In order to understand better the problem of the boundary T = 0, it is useful to recall the equivalence between the 
equation w = and the so-called Mayer-Lie system of partial differential equations [herein, stays for any extensive 
variable different from U and ^, for the corresponding intensive variable] 

%{X\...,X"+^)=^i{U,X\...,X^+') fori=l,...,n+l. (64) 



dXi 

One can also assign an initial condition 

U{Xl...,XS+') = Uo (65) 

and thus define a Cauchy problem for the above Maycr-Lic system. The integrability condition u A doj = in the 
inner part of the manifold is sufficient for a Pfaffian form in order to ensure the existence and the uniqueness of 
the above Cauchy problem. This means that the Cauchy problem with initial point on the T = boundary allows 
solutions which lie in T = 0. If w is also on the boundary, then it is evident that the aforementioned curves are the 
only possible solutions to the above Cauchy problem with initial point on the surface T = 0. In other terms, if a; G 
everywhere, then T = is a leaf of the thermodynamic foliation. But, a priori, one can consider also a PfafHan form 
CO such that it is continuous on the boundary T = but non-necessarily there. The uniqueness of the solution 
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of (64) with initial condition on the surface T = could be ensured if the functions £_i{U,X^, . . . ,X"~^^) are locally 
Lipschitzian with respect to U uniformly with respect to , . . . , X"^^ in a neig hborhood of {Uq, X^,..., X^+^). If 
even this condition fails, then the continuity of ui, i.e., the continuity of ^ also in T = can allow multiple solutions 
of the differential equation (64). 

Let us consider the following differential equation which describes isoentropic curves in the special coordinated 
adapted to the boundary introduced in the previous subsection: 

^ = Y: UB{r),X\r), . . . ,X"+i(r))^ (66) 

i 

This equation can be easily obtained from (64). Let us consider at least piecewise functions X^{t), . . . , X"+^(r) 
for T G [0, 1]. These hmctions are arbitrarily assigned. Let us consider a solution curve such that lim,-^T-o B{t) = 
for To € [0, 1]. By hypothesis, X-^{t), . . . ,X"+^(t) are finite for t — >■ tq. Then, by continuity, such a solution can be 
extended to r = tq, i.e. T = cannot be a leaf. This happens as a consequence of well-known theorems on the ordinary 
differential equations, sec [45], pp. 67-68. Let us then consider a point (B(0) = 0,^1(0) = X^, . . . ,X"+i(0) = X^+^) 
on the surface T = 0. In each neighborhood of this point, one can find inner points, each of which belongs to a surface 
S = const. In fact, each point of the surface B = is a limit point for the nearby inner points of the thermodynamic 
domain and each inner point has to belong to a S = const, integral manifold, because of the integrability condition. 
If T = is a leaf [or if the connected components of T^^(O) are leaves if T = is not connected], the only possibility 
is that, in approaching T = 0, one is forced to change leaf S — const., i.e., it is not possible to approach T = by 
remaining on the same leaf S = const., otherwise the inner solution of w = could be extended to T = 0. Notice 
that, in case there exist two solutions of the Cauchy problem for the differential equation (66) with initial condition 
(5(0) = 0,^1(0) = Xo\...,X"+i(0) = XJ+^), one lying in the T = surface and the other leaving the T = 
surface, then these solutions are tangent at the initial point. This means that, in case of existence of the limit as 
S — > 0+ of the entropy, when (N) is violated, there are surfaces S = const, which are tangent to the submanifold 
T = 0. We now prove these statements. 



1. validity of (N) 

If (N) holds, then T = plays at most the role of asymptotic manifold for the inner leaves S = const. No inner 

integral manifold can intersect T = 0, i.e., T = is a leaf. In order to approach T = at finite deformation 
parameters, one has necessarily to change from one adiabatic surface S — const, to another one, it is impossible to 
approach T = by means of a single adiabatic transformation. We have then a mathematical explanation of the naive 
unattainability picture sketched by means oi S — T diagrams one finds in standard textbooks on thermodynamics (cf. 
also the definition of P{x) transformations in [2]). 



2. violation of (N) 

If, instead, (N) is violated and lims^o+ S{B, Xq, . . . , X^'^^) = 5(0, Xq, . . . , Xq +^), then the inner integral manifold 
< 5(0, Xl, . . . , Xq'^^) = const, exists and can be continuously extended to T = 0. 
This can be proved by means of a variant of the implicit function theorem. For simplicity, we put here 

X = XS...,X"+^ (67) 

Let us consider a point (0,Xo) which is not a local minimum for S, i.e. it is not such that 5(0, Xq) < S{B,X) in 
a neighborhood of (0,Xo). Such a point surely exists if (N) is violated, as it is easy to show (cf. also [19]). We are 
interested in the zeroes of the function 

a{B, X) = S{B, X) - 5(0, Xq). (68) 

The fimction (t{B, X) is a continuous function which is monotonically strictly increasing in B everywhere in the 
domain V U dV, because S{B,X) is, by construction, a strictly increasing monotone function in B, as it is evident 
from (54). 

In particular, we wish to know if there is a continuous function B(X) defined in a neighborhood of {0,Xq) such that 
a{B{X),X) = and such that B{Xo) = 0. If (0, Xq) is a strict local minimum for 5, then by definition there exists a 
neighborhood of (0, Xq) where cr{B, X) > 0, thus the aforementioned function B{X) does not exists. If it is a weak local 
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minimum, in the sense that S{0, Xq) < S{B, X) in a neighborhood and the cquahty is aUowed, again a{B{X), X) = 
does not admit sohitions. in fact if S{B, X) — S{0, Xa) is allowed in a neighborhood W D [0, Bq) xVoi (0, Xq), with 
V B Xq open set, then S{0,X) < S{0,Xo) because S{0,X) < S{B,X). Thus, being (0, Xq) a local minimum, one 
can find a smaller neighborhood where S{B, X) = S{0, Xq) is impossible for any B > {S{0, X) = S{0, Xq) is instead 
allowed). 

Notice that (0, Xq) cannot be a local but non global minimum for S under the natural requirement that each surface 
S = const, corresponds to a unique integral manifold of SQrev [this means that isoentropic states are path-connected. 
Cf. [39]]. In fact, in homogeneous thermodynamics, there exists an integral manifold of SQrev such that S{B,X) = 
S{0,X()) ~ const, if {0,Xq) is not a global minimum, and to such an integral manifold (0, Xq) would not belong if 
(0,Xo) is a local minimum. Notice also that no point belonging to the boundary B = can be a local maximum, 
because S{B, X) is a strictly increasing monotone function in B. 

In any convex neighborhood W of (0,Xo) there exist (_B+,X+) and {B~,X~) such that a{B~^,X^) > 
and (7{B^ , X^) < 0, because (0, Xq) is not a local minimum. By continuity, for any convex neighborhood of 
(0, Xo) there exists such that a{B°,X") = 0. The point (0,Xo) is then a limit point for the set 

Z{a) = {{B'^ , X'^)\a{B'^ , X^) = 0}, which is a closed set because a is continuous. In order to show that a solu- 
tion B{X) for a{B{X),X) =0 exists and is unique, we introduce the following auxiliary function 



a{B,X) 



a{B,X) forS>0; 
5(0, X) - S{0, Xo) + B for B < 0. 



This function a{B,X) is a continuous function which is monotone strictly increasing in B also for B < 0. We have 
extended then a to negative values of B, which is shown to be an useful trick. We cannot use the standard form of 
the implicit function theorem because dS/dB diverges at B = 0. Nevertheless, the proof is a variant of the standard 
proof of the implicit function theorem for scalar functions. We have that (t(0, Xq) = 0. By the monotonicity property 
one has that there exist Bi < < B2 such that ^{Bi.Xq) < < (7{B2, Xq). By continuity, there exists an open 
neighborhood Rb Xq such that a{Bi,X) < < a{B2, X) for all X G R. Then, from the intermediate value theorem 
it follows that there exists a value B G {Bi^Bi) such that a{B,X) = for any fixed X G R. Monotonicity ensures 
that B is unique for each fixed X E R. The function B{X) : i? ^ IR is then defined as the map defined by B{X) = B, 
where B is the solution of a{B,X) = for each fixed X G R. Such a function satisfies B{Xo) = and is also 
continuous in X. This follows again from the fact that a is a continuous function which is monotonically strictly 
increasing in B. Being continuous, one has that the set 

Z{a) = {{B,X)\a{B,X)=0} (70) 

is a closed set. Given a sequence {X„} c R such that X„ X G R iov n ^ 00, one finds that there exists a 

unique (by monotonicity) B such that {B,X) G Z{a). This means that lim„^oo B{Xn) = B{X) = B, i.e., B{X) is 

continuous. 

It is evident that 

Z{a) 2 Z{a) (71) 

and we have to get rid of the spurious solution B = —{S{Q,X) — 5(0, Xq)) < which could occur for 5(0, X) — 
5(0, Xq) > 0. But this solution cannot hold for any X G R, because in any neighborhood of (0, Xq) there exist points 
(0, X~) such that 5(0, X~) — 5(0, Xo) < 0. Thus, being the spurious solution not defined in the whole R, the above 
theorem allows to conclude that B{X) > surely exists. 

Notice that, for inner points B(X) > 0, this solution is actually a leaf of the foliation defined by the integrable at 
least Pfaffian form lo. As a consequence, for inner points B{X) is at least C^. One can also calculate the gradient 
of B(Xi,...,X"+i): 



dB dB 



= (72) 



^axi''"' ax"+\ 

where the latter equality is due to the fact that B{X) satisfies the Mayer-Lie system. It is then evident that 

dB dB \ 

axT'--'ax^j^o 

i.e. B{X) reaching B = is tangent to iJ = 0. 
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The above condition about the absence of inner integral manifolds arbitrarily approaching T = is also sufficient. 
Notice also that it does not forbid the inner leaves to asymptotically approach T = as some deformation variable, 
say X*, is allowed to diverge: \X'^\ — > 00 as B — > 0+. The unattainability is clearly ensured because of such a 
divergence. Let us consider the following equation: 

^ Ub, X\B), . . .,X-+\B)) ^-1 = 0, (74) 

i=l 

which is another rewriting of the above equation where B plays the role of independent variable and where n functions 
X*, say X^{B)^ . . . ,X"(_B), are arbitrarily assigned. For simplicity, let us put X'^'^^[B) = X[B). One could also 
consider X^ = Xq, . . . , X" = Xq = const. Then, one obtains 



dX 1 



(75) 



It is evident from our discussion above that any solution X{B) has to be such that |-'i^(iJ)| 00 when iJ ^ 0+ if 
T = has to be a leaf. Some examples are given below. 



3. examples 

Let us consider 

u = dU+^^dV; (76) 

the domain is chosen to be < ?7, <V and the PfafRan form ui is everywhere. One has 

f=lu (77) 

which vanishes for U = 0. The boundary U = is an integral submanifold of w. Let us consider the Cauchy problem 

^ = -^^ (78) 
dV 3V ^ ' 

UiVo) = 0. (79) 

It is evident that the only solution of this problem is U = 0, which is a leaf of the thermodynamic foliation. By 
integrating, one finds the (concave) entropy S = cq U'^I^V'^I^ and T = 5/(3co) (U/V)"^^^ [cq is an undetermined 
constant]. (N) is satisfied. Along an isoentropic surface So > 0, one finds 

U=(^^y^\-V^ (80) 

and T = 0, i.e. U = can be approached only for V ^ 00. 

Let us consider a PfafRan form having the same domain <U, <V 

/f^x 2/3 

iv = dU+l-] dV; (81) 

the Pfaffian form w is not on the boundary {7 = 0. One has 

f = U + ;72/3y 1/3 (82) 
which vanishes for U = 0. The Cauchy problem 



dU fU\ 



2/3 



dv = -[v) 

U{Vo) = (84) 
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allows two solutions: 

U = (85) 

and 

U = {Vo^/^ - V^/^). (86) 

The latter solution holds for < ^ < Vb, and it can be easily identified with the isoentrope S — Sq = cqVo, where 
S = Co(J7^/^ + V^^^)^ is the (concave) entropy. (N) is violated and the two solutions are tangent for U — 0. 

Let us consider the following example, which is inspired to the low-temperature behavior of a Fermi gas. The 
PfafRan form one takes into account is 



where c is a positive constant. This PfafJian form is integrable and the integrating factor is 



Ar5/3 

= 2U -2c 

Then the zero of the integrating factor occurs for 



/ = 2t^-2c^- (88) 



Ar5/3 

U = c^, (89) 

and, by construction, being / > 0, one imposes U > b{V,N) = c N^/^/V"^^^. The function b{V,N) is extensive and 
convex. Let us define 

B = U-c^; (90) 
This coordinate transformation is regular in S = 0. We have 

p{B,V,N) = ^^ (91) 
jl{B,V,N) = -^^, (92) 

and 

f = 2B. (93) 

Then, one finds 

S = a BV2yi/37vi/6^ (94) 
(a is a proportionality constant) which can be easily re-expressed in terms of {U, V, N): 

S = a (c/F2/3ArV3 _ civ2) . (95) 
Notice that, along S = So = const, one has 

which can approach B = only for F — > oo and/or N ^ oo. Moreover, 

T=l (c/y2/3^1/3 _ ^^2^ ^-1/3 y-2/3 (97) 

and dT/dU diverges as T ^ 0+. 
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H. conditions for the validity of (N) 



We can show, by means of purely thermodynamic considerations, that: 

in Gibbsian variables, if the homogeneous Pfaffian form is (at least) also on the boundary Z{f) and if the entropy 
S is concave, then (N) holds. 

{sufEcient but not necessary condition) 

Recall that S' > by construction and also on statistical mechanical grounds. Let us assume that the Pfaffian 
form SQrev is everywhere, also on the boundary, and that it satisfies the integrability condition in the inner part 
of the manifold. Then, the integrating factor / is a function everywhere. 
In variables U, V, N one has 



where p, fi are . Consider 



f = U+ p{U, V,N)V- ii{U, V, N) N (98) 



df dT 



S cannot be non- negative and concave near T = if it diverges as T — > 0+ [2]. See also sect. VII. Then concavity 

and positivity also near T = force the entropy to be finite in the limit T ^ 0+. Moreover, one has dT/dU — 
1/Cjfi_ ^jfn+i, where Cx^ ,...,x^+^ is the standard heat capacity at constant deformation parameters, which has to 
vanish in the limit as T ^ 0+ because S has to be finite in that limit. As a consequence, dT/dU = ^ oo 

as T ^ 0+. Then, by inspection of (99), it is evident that / is also on the boundary only for S* ^ 0+ as T ^ 0+. 
As a consequence of this theorem, we can conclude that any violation of (N) is involved with a Pfaffian form that is 
not also on the boundary, as it can be easily verified by considering the examples violating (N) in section VII. 



1. geometrical aspects of f = 

If SQrev is at least of class everywhere, then d{SQrev) is continuous and finite on the surface / = (that is, 
T = 0). From a geometrical point of view, the surface / = is a so-called separatrix according to the definition of Ref. 
[40], in the sense that the 2— form {SQrev) A df vanishes with / on this surface. One has (SQrev) ^df = —f d{SQrev), 
as it can be easily verified by using (21) and standard Cartan identities. / = is also a surface where the symmetry 
associated with the vector field Y becomes trivial [41], in the sense that it becomes tangent to the submanifold / = 0. 
In the case where the complete integrability of SQrev is preserved also on the surface / = 0, one has that / = 
appears as a special leaf of the thermodynamic foliation, it is indeed the only leaf which is left invariant by the action 

of y. 

In general, for a sufficiently regular w at / = 0, the property for / = to be a separatrix in the sense of Ref. [40] 
ensures that the solutions of SQrev = in the inner part of the thermodynamic manifold cannot intersect the integral 
submanifold / = 0. Then, the integral manifold / = is a leaf of the foliation defined by SQrev if it is also a separatrix. 
Moreover, it is a special leaf, being a separatrix. 

We have shown above that 5 0"*" as / ^ 0"*" for a; e everywhere. We can note that, if the condition for SQrev 
to be also on / = is relaxed, then the condition 

\ lim^ f d{SQrev)\ = (100) 
is not sufficient for ensuring (N). A simple counterexample is given by 

u) = dU+ (^^y dV (101) 
where <U, < Vi < F and where / du vanishes as / — >■ 0+ but (N) is violated. 

The requirement that w is everywhere is, at the same time, too restrictive. In fact, it is evident that the following 
homogeneous Pfaffian form, defined for D = {0 < J7 < V/e^, V > 0}) 

oj = dU + p{U,V) dV, (102) 
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where 

r ^(-l-log(^)) iiQ<U<V/e^ 
p{U,V) = { 

[ if [/ = 

satisfies (N) but is not of class at T = (the corresponding entropy is S = —V/ {\og{U /V)), which is concave on T) 
and can be continuously defined to be zero when [7 = 0). Moreover, this example shows that even the more general 
setting 

S concave, ^ finite as T ^ 0+ =^> 5 ^ 0+ as T ^ 0+ 

(which is trivial, because d f /d U = 1 + S {d T/d U)) does not correspond to a necessary condition for (N). 



2. condition (HOM) 

In order to give a necessary and sufficient condition for (N), we use the following interesting property of dQrev We 
have 

SQrev _ dS nr\-i\ 

Let us consider J SQrev/.f, where 7 is a curve having final point at temperature T. If (N) holds, then 
SQrev/.f —00 as T ^ 0+. In fact, if (N) holds, whichever the path 7 one chooses, the integral of dS/S 
diverges to —00 as T ^ 0+. 

If, instead, SQrev/f —00 as T — > 0"*" whichever path is chosen for approaching T = 0, then 5 — > 0+ in the same 
limit. 

Then the following theorem holds: 

(N) J SQrev/ f — * —00 as T ^ 0+ whichever path is chosen in approaching T = 0. 



{condition (HOM) in the following) 

Notice that, because of the concavity of S, one cannot have J^^ SQrev/.f +00 as T ^ 0+, because a non-negative 
and concave entropy cannot diverge (cf. sect. VII). Thus, once ensured the concavity property for S (cf. [39]), one 
has only to check if the above integral diverges along any rectifiable curve approaching the surface T = 0. 

Notice that, in this form, the above theorem allows to neglect the problem of the actual presence of the boundary 
T = in the physical manifold. This formulation is also coherent with the fact that (N) is formulated as a limit for 

If (N) is violated and the limit limT_,o+ S exists, then (103) is integrablc^ along any path approaching T = with 
positive entropy (it is not integrable along any path approaching T — with vanishing entropy). For example, let us 
consider the following toy-model Pfaffian form 

- = ^^+U+"°(yl^) '^+U+^°(FAOV^) ^'''^ 
where ao > 0, /3o < are constants and the domain is restricted by V/N > —Po/cxo- Then one has 

The integrating factor / vanishes as C/ ^0. oj/f is integrable along any path such that V/N > — /3o/ao, in fact, if 
g{V,N) stays for a positive function, one has / ~ C/^/^ g{V,N) as U ^ 0+. If V/N = -/3o/ao, then / ~ 3 J7 as 
?7 — > 0+ and w/f is no more integrable near U = 0. Notice that the entropy which corresponds to this Pfaffian form 
is S = 3{UV N)^^^ + aoV + I3a N >0. (N) is violated and S vanishes on the submanifold U = 0, V/N = -Po/cto- 
Notice that, if S is not concave but simply positive, then condition (HOM) is still equivalent to (N). 
Notice that (HOM) is not affected by the connectedness properties of / = 0. 
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I. a sufficient condition for a continuous S at T = 



The entropy S is required to be continuous at T = 0. This can be obtained as follows. It is evident that the 
continuity of S implies the continuity of S. Then, one can impose conditions which allow S to be continuous at T = 0. 
Nevertheless, the continuity of S excludes, by direct inspection, the possibility to obtain 5 = at T = 0. Thus, 
one has to find further conditions on S in order to allow the possibility to get a vanishing entropy at T = 0. Let us 
consider 

/■^ . 1 

S(B,X\...,X''+^)= / dB — . +S(Bo,X\...,X''+'^). (106) 

We wish to know if the limit as S — » 0+ of 5 exists. A sufficient condition is the following: there exists a positive 
function <p{B) such that 

< 0(B) V B e (0, Bo] and . . . , e C, (107) 



/(B,Xi,...,X"+i) 

where C C IR"+^ n V is any open bounded set contained in V and 

lim / dB 4>{B) < 00. (108) 

Then J?^ dB 1// is uniformly convergent and, being S continuous for B > 0, one finds that S can be extended 

continuously also at B = 0. This condition ensures that (N) is violated. The above condition does not leave room for 
S" = for some (but not all) values of X^, . . . , e V. Actually, continuity on a open bounded subset R C IR"'''^n'E> 
of the allowed values for the variables X^, . . . can be also obtained by assuming that (107) holds on R and 

not for any open bounded set C contained in V. In this case, S is continuous at T = for X^, . . . , X""*"^ e R. 
In order to obtain a condition ensuring (N) a sufficient condition is the following: there exists a positive function 
^{B) such that 

— . ^- >(I)(B) VBe(0,^o] and . . . , e C, (109) 

f{B,X\...,X^+^) ^ ' ^ ^ ^ ' 

where again is any open bounded set, and 

/ dB (l){B) = -00. (110) 



lim 

B^0+ J Bo 



Then, because the above divergence of the integral is uniform in X^ , . . . , one finds that S = log(S') —oo as 

T ^ 0+, i.e. 5 — > 0+ as T — > 0+. Even in this case, one can allow the function ^(-B) to be different for different 
subsets whose union covers all the values of the deformation parameters. 



1. (N): integral criterion 

We can consider 



^ = S{B, X\..., X"+^) - SiBo, X\..., (111) 
r / 

due to the singularity in the integrand in / = 0, when one considers a path approaching T = the integral has to be 
intended as improper integral. Nevertheless, according to a common use, we bypass this specification in the following. 
We can write 

S(B,X\...,X''+^)-S(Bo,X\...,X^+^)= / dB } : (112) 
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We can sec that, if w e C^{V U dV), then the above integral diverges in such a way that 5 ^ 0+ for B — > 0+. In 
fact, one has 

i,{B,X\..., = ^k{B,X\..., - a(0, X\ . . . , = ^{0,X\..., B + 0{B^) (113) 

oB 

and 

/ = fc(0,X\...,X"+i) 5 + 0(^2), (114) 
where surely fc(0, X^, . . . , X"+^) > because / > 0. As a consequence, near B = the integrand behaves as follows: 

UJ 1 

-adB. (115) 

f kB 

Then the integral diverges as log(S) for B — > 0+. The function (j>{B) is 

<^(B)= sup (^(0,Xi,...,X"+i)) 1, (116) 

where C is an open bounded set contained in D. 



J. inaccessibility (C) and the failure of (N) 



If (N) holds, the surface T = is adiabatically inaccessible along any adiabatic reversible transformation starting 
at T > 0, and it is a leaf of a foliation. There is no isocntropic surface reaching T = and the property (U) of 
unattainability is automatically ensured, as well as the principle of adiabatic inaccessibility (C). 

The violation of (N) is instead very problematic from the point of view of (C) and of the foliation of the thermody- 
namic manifold. If (N) is violated, T = is not a leaf and it is possible to reach T = along inner (would-be) leaves 
S = const. Actually, one does not find a foliation of the whole thermodynamic domain; if T = is included in the 
thermodynamic manifold, one finds an "almost-foliation" , i.e. a foliation except for a zcro-mcasurc manifold, in the 
sense that to the proper inner foliation generated at T > is joined a integral manifold T = (the adiabatic boundary 
of the thermodynamic domain) which breaks the adiabatic inaccessibility, even if only along special paths passing 
through T = 0. In the spirit of the thermodynamic formalism, we agree with Einstein's statement that the existence of 
such adiabatic paths is "very hurtful to one's physical sensibilities" [23]. It is also evident that the Carnot-Ncrnst cycle 
discussed in sect. IV A is allowed, unless some discontinuity occurs or the thermodynamic formalism fails according 
to Planck's objection, and that the objections against its actual performability can hold only in restricted operative 
conditions (from a mathematical point of view, a path contained in the surface T = is different from a isoentropic 
path at T > reaching the absolute zero of the temperature). Moreover, the approach to the problem by means of 
SQrev reveals in a straightforward way aspects which other approaches cannot easily point out. 

A further remark is to some extent suggested by black hole thermodynamics, where (N) is violated but states at 
T = have S = 0. Cf. [20] for a study in terms of Pfafhan forms. In order to avoid problems occurring with the 
surface T = if (N) is violated, one could introduce a further hypothesis. One could impose that the entropy is 
discontinuous at T = 0, and that 

So < inf S(F,X\...,X"). (117) 

v,x\...,x^ 

One could then impose that Sq = for all the systems, which would allow to recover an universal behavior. Even 
the adiabatic inaccessibility would be restored, because the second law would inhibit to reach T = adiabatically. 
This behavior characterizes black hole thermodynamics. Concavity would be preserved, as well as superadditivity. 
However this choice is arbitrary and even unsatisfactory, because a well-behaved foliation of the thermodynamic 
manifold is obtained by hand by means of the discontinuous entropy S just constructed. In fact, the foliation of the 
thermodynamic manifold, if (N) is violated, is obtained as the union of the usual foliation at T > and a special leaf 
at T = 0. This foliation is generated by a Pfaflian form only in the inner part of the manifold. 
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VII. NOTES ON GIBBSIAN APPROACH 



We recall that in Gibbsian approach [47] , the existence of the entropy is a postulate, because the entropy appears 
in an axiomatic framework. Sec also Rcfs. [48,49,17]. In a certain sense, very loosely speaking, Gibbs starts where 
Caratheodory leaves [50]. This can be considered the reason why in Gibbsian approach the problems which can be 
associated with the surface T = as in the previous section appear to be less evident. 

Let us assume the Gibbsian approach to thermodynamics, and write the so-called fundamental equation in the entropy 

representation: 

S = S{U,X^,---,X''+^) (118) 

where , ■ ■ ■ . X"'~^^ arc extensive deformation variables and U is the internal energy. S is required to be a first order 
positively homogeneous function and, moreover, a concave function (for mathematical details about convexity we refer 
to Refs. [51,52]). The former property ensures the extensivity of the entropy, the latter ensures the thermodynamic 
stability of the system against thermodynamic fluctuations. The second law AS > for an insulated system is also 
ensured. 



A. extension of 5 to T = 

Let us define 

I{U, X\---, X"+i) = -S{U, X\---, X"+i). (119) 

In what follows, x stays for a state in the thermodynamic manifold: x = {U, X^, • • • , X"+^). The function I(x) is, by 
definition, a convex function and positively homogeneous function. As a consequence, its epigraph is a convex cone. 
This convex function I is defined on C [we change symbol for the domain, what follows holds for a generic convex 
function in a generic convex domain]. There is a preliminary problem. One has to define I{x) on the boundary dC 
and obtain again a convex function. This is made as follows [52] : / can be extended to the set 

= C U dCf, 
where 

dCf = {yedC\ liminf /(a;) < oo} 

and 

/(y) = liminf /(x) WyedCf. (120) 

x^y 

The above extension is convex on a convex set. In general, one cannot substitute lim inf /(.t) with lmix^yl{x) 
because the latter may not exist [52]. Moreover, the behavior of the convex function / = —S at the boundary has to 
be such that 

liminf I{x) > -oo (121) 

for any xq belonging to the boundary of the convex domain [cf. problem F p. 95 of Ref. [52]]. 

In the case of S, then the non-existence of the above limit can be considered unphysical. In fact, it can also mean 

that the entropy could approach a different value for the same state along different paths starting at the same initial 
point. In the latter case, its nature of state function would be jeopardized, it requires at least the existence of the 
limit, that is, the independence of the limit from the path chosen. On this topic, see in particular [2,3]. In any case, 
the definition offered by the theory of convex functions 

S{y) = lim sup S{x) V y e dCf (122) 

x^y 

is a rigorous formal prescription, but it is not clear to the present author if it could be relevant to the physics at hand, 
if the limit does not exists. 

Then, we assume that S admits a limit for each point of the boundary T = 0, thus 

S{y) = lim S{x) V y e {T = 0}. (123) 
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Under this hypothesis, we can extend uniquely S at T = 0. The surface T = represents (a part of) the boundary 
dT> for the domain, then it belongs to the closure of the convex open set V. A convex set is dense in its closure. As 
a consequence, a continuous function G defined in V can be uniquely extended by continuity on the boundary dV if 
(and only if), for each point Xb € dV the limit 

lim G 

x&'D—^xt, 

exists. 

It is still to be stressed that, for a non-negative concave S, one has to find \imx^o+ S < oc as a consequence of (121). 

B. attainment of the lower bound of S 

Gibbsian approach allows us to conclude immediately that, if the upper bound Iq of / is attained, then it has to 
be attained on the boundary of the thermodynamic manifold (if a convex function / should get a maximum value Iq 
in a inner point of its convex domain, it would be actually a constant function in its domain) [52] . Then, if the lower 
bound So of S is attained, it is attained on the boundary of the domain of S. Moreover, under very simple hypotheses 
on the domain, the upper boimd of / is actually attained [51,52]. In particular, it can be attained at an extreme point 
of the boundary. We recall that an extreme point of a convex set is a point belonging to the boundary of the set such 
that it is not an inner point of any line segment contained in the set. For example, if the set is a closed rectangle, the 
extreme points are the four vertices; if the set is a circle, all the points of the boundary (circumference) are extreme 
points. But notice also that, being the thermodynamic domain a convex cone, there is no extreme point apart from 
the origin of the cone (which cannot be considered a physically meaningful state [39]). 
Note that, given the surface T = 0: 

r([/,X\...,X"+i) = 
U = Uo{X\...,X^+^) 

then, for each point on this surface, as a consequence of the homogeneity of the entropy, it holds 

Six Uq{X\ X"+i), A X\...,X X"+^) = A SiUoiX\ X"+^),X\ . . . ,X"+i) 

At the same time, the intensivity (i.e., homogeneity of degree zero) of T implies 

T(A Uo{X\ X"+i), XX\...,X X"+i) = T{Uo{X\ X"+i), X\ . . . , X"+i) = 0. (127) 

Then, if Uo{X'^, • • • , X"+i), X\ . . . , = Xg,a = 0, . . . ,n ave the points belonging to the surface T = 0, the cone 

Ko = {XS I A XS €Ko, X> 0} (128) 

is contained in the surface T — because of the intensivity of T. As a consequence, in case of violation of (N), one 
could find a system at T = having an arbitrarily high entropy. Only if S = at T = this cannot happen, because 
5 = is a fixed point under scaling of the entropy. 

C. values of S at T = and the hypothesis of multi-branching 

We have assumed a continuous S' at T = 0. From the point of view of the Landsberg's discussion about a multi- 
branching near T = 0, we have then simply to discuss the following topological problem. Is the set Z{T) a connected 
set? In the case it is connected, then we can surely conclude that no multi-branching can occur near T = 0. In fact, 
the range of a continuous function on a connected set is a connected set, that is, the range of S* at T = is a connected 
set contained in IR. It has to be an interval (violation of (N)) or a single point (validity of (N)). It is interesting to 
underline that, even if the set Z{T) is not connected, (HOM) ensures the validity of (N) and viceversa (there is no 
possibility to find two branches like the ones in Fig. (la), because both have to start at 5 = 0,T = 0). A possibility 
for getting a multi-branching is to violate the concavity at least near T = 0. For an interesting example see Ref. [20]. 
[Another possibility to get a multi-branching could be to consider a system allowing for states at T < 0, but in this 
case two distinct branches would be found on two difi^erent sides of T = 0] . 
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(124) 



(125) 



(126) 



D. violation of (N) and Landsberg's hypothesis 



Consider the following toy-model: 

S = jo u^+SoV (129) 



where 70 > 0, (Jq > and < a < 1. Then 



U 



Jo a \V 



(130) 



which vanishes as U ^ 0+: T = <^ ?7 = 0. For the domain let us consider !F = {U > 0} U {V > Vq}. I is 
maximum, that is, S is minimum, at the extreme point (0, Vq), as it is evident. [Notice that in this example the 
domain is not a convex cone because we introduce a lower bound Vq for V, as it is physically reasonable in order to 
justify thermodynamics on a statistical mechanical ground. If one considers = {U >0}L}{V > 0}, then inf (5) = 0, 
which is approached at the only extreme point (0,0) of the cone]. If 5 > 0, then (N) is violated and S can assume 
an interval of values at T = 0. If So = 0, then (N) is satisfied and / is maximal on the line U = 0. A special case is 
represented by the photon gas, where a = 3/4. 

This toy-model corresponds to the following behavior of 5 as a function of T and V: 

S{T, V) = (eo TT^ + So) V (131) 
which, for Sq ^ 0, violates (N). It is useful to pass to the energy representation 

The domain is Q = {Vq < V < S/So}- The T = surface corresponds to V = S/Sq. We have 

T=-^ {S - So V)^-' ; (133) 



the pressure is 



p=- (—] " {S-So y) V-i {S{l-a) + Va Sq). (134) 
a V70/ 



The isoentropic S = So has equation 

U^y^^l^So^^y (135) 

and reaches T = when V = Sq/So during an adiabatic expansion. It is easy to see that it is tangent to the T = 
surface. The adiabatic expansion has to stop there, because of the structure of the domain. One can wonder if 
any physical reason for such a stopping exists. It is useful to come back to Landsberg's suggestion about a possible 
vanishing of the adiabatic compressibility: 

In our case, we get 

Ks = a^4 ^^V^ (S-SoV)"^ (137) 

and three cases occur: when 1/2 < a < 1 then Ks ^ as T ^ 0+, in such a way that the elastic constants of 

the system diverge in that limit, forbidding any further expansion (Landsberg's behavior); when < a < 1/2 then 
Ks ^ 00 as T ^ O"*", the elastic constants vanish and the behavior of the system is pathologic (the system appears to 
be "totally deformable" in that limit); when a = 1/2 then Kg — > 7o V'^/{2 S"^) which is in any case finite {S is surely 
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positive and not zero) and a physical hindrance against reaching T = is not apparent. For a = 3/4 the behavior for 
5 as T ^ 0"*" is hke the one of Debye model, except for the shift SqV which allows the violation of (N); the adiabatic 
expansion stops at T = because of a vanishing Ks- For a — 1/2 one obtains a behavior similar to the one of an 
electron gas near T = 0, except for the shift SqV; no vanishing or divergence of Ks is allowed at T = 0. 
We wish to underline that, if one consider the entropic fundamental equation of an electron gas at low temperature 
and shifts it by Sq V, then the situation is still different, because of the zero-point mode contribution. In fact, such a 
contribution allows to get a positive and finite p as T ^ 0+ and, moreover, a positive and finite Ks,n as T — » 0+ (the 
domain is 1^ > 0, A/' > 0, f7 > Uo{V, N), where N is the particle number). As a consequence, Landsberg's mechanism 
does not seem to be available. Instead, even if the zero-point mode contribution is taken into account in the case of 
a Debye crystal, the pressure is still positive as T ^ O"*" but Ks,n vanishes. 

E. further properties 

Another point that can be underlined is the following. Let us assume to extend the convex function / = — S' to 
all of IR", by defining / = -|-c» outside its domain dom /. Then, replace this function with its closure. This is the 
same procedure which is prescribed in Ref. [53] for the internal energy U. Then 7 is a closed proper convex function 
which is essentially smooth, that is, jV J| ^ oo for any subsequence converging to a boundary point. In fact, in the 
gradient of / appears the factor 1/T which diverges as the boundary T = is approached. This allows to obtain in 
thermodynamics a convex function of Legendre type, which is relevant for the discussion of Legendre transformations 
in thermodynamics. The difference of the entropy representation with respect to the energy representation is evident 
from this point of view. In fact, the opposite conclusion appears in Ref. [53] for the energy representation: no Leg- 
endre type hmction exists in thermodynamics. But we think, from the discussion of the previous section, that the 
entropy representation is more fundamental with respect to the energy representation at least for the discussion of 
the boundary T = 0. 

Summarizing our analysis in Gibbs framework: 

gl) if > , then S = can be attained at a point of the boundary of the domain; 

g2) it could be that S* > at other points on the surface T ~ 0. Then, homogeneity implies that, by scaling, a system 
with an arbitrarily high zero temperature entropy could be obtained; 

g3) models exist where the violation of (N) does not imply the attainability of T = and the violation of AS' = for 
adiabatic reversible transformations (states at T = are not available and so no such violation can occur at T = 0). 
But these systems display a non- universal behavior (i.e., a behavior which does not appear in other models). 

VIII. CONCLUSIONS 

We have discussed the status of third law of thermodynamics and we have given an heuristic argument in favor of 
the entropic version of the third law. Then, we have analyzed the law both in Caratheodory's approach and in Gibbs' 
approach to thermodynamics. 

In particular, Garatheodory's approach shows that for T > the thermodynamic manifold can be foliated into 
leaves which correspond to isoentropic surfaces. The only hypothesis is that the Pfaffian form SQrev is integrable and 
in the inner part [T > 0) of the thermodynamic manifold. At the boundary T = 0, which is assumed to be an 
integral manifold of the Pfafhan form SQrev , the aforementioned Pfaffian form is allowed to be also only continuous. 
The special integral manifold T = is problematic from a physical point of view, because it can also be intersected 
by the inner (would-be) leaves S = const. In the latter case, (N) is violated and one obtains an almost-foliation of 
the thermodynamic manifold, where the inacx'cssibility property fails, even if only along special adiabatic paths which 
pass through the surface T = 0. For an entropy which is continuous also at T = 0, (N) holds if and only if T = is a 
leaf. This is a remarkable result, the validity of (N) is strongly related to the possibility to obtain a foliation for the 
whole thermodynamic manifold, including T = 0. We have shown that, if the Pfaffian form is everywhere, then 
(N) is preserved. Physical assumptions and mathematical conditions have been discussed. 

In another paper [19] , further conditions leading to the third law are discussed. 

We add herein some notes about the conditions ensuring (N) in quasi-homogeneous thermodynamics introduced in 
[54]. Also in the quasi-homogeneous case (N) holds iff lim7-^g+ S = 0. The analysis of sects. VI F and VI G hold 
with obvious changes; moreover, condition (HOM) holds unaltered, and, if the Pfaffian form oj is everywhere, then 
(N) holds (this can be shown by using a criterion analogous to the one appearing in (VII 1). 
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